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We look for a Brans-Dicke type of generalization of Horava-Lifshitz gravity. It is shown that 
such a generalization is possible within the detailed balance condition. The resulting theory reduces 
in the IR limit to the usual Brans-Dicke theory with a negative cosmological constant for certain 
values of parameters. We then consider homogeneous and isotropic cosmological situation in the 
context of this generalized theory, and find some interesting features of the Brans-Dicke scalar field 
in determining the behavior of the universe. 
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Recently, a new theory of gravity has been proposed by Horavafll-Q. This theory, being based on anisotropic scaling 
of space and time, breaks the spacetime symmetry. It has a much better UV behavior than the theories with spacetime 
diffeomorphism symmetry such as general relativity, but reduces to Einstein's gravity in the infrared limit, thereby 
recovering spacetime diffeomorphism symmetry. Physical constants such as the speed of light, Newton's constant, 
and cosmological constants all emerge from the relevant deformation of the non-relativistic theory at short distance. 
These interesting features as well as related findings have received a great deal of attention 

On the other hand, even if we consider only the low energy limit of the gravity there are many alternative theories 
and extensions of the Einstein theory. In particular, in the context of cosmology various models with a scalar field 
have been considered and a possible role of the scalar field in explaining the behavior of the universe in the early 
inflationary stage as well as the late stage has been investigated [3, @|. Therefore, it would be interesting to see if 
Horava's theory can be extended in such a way that in the infrared limit it reduces to those alternative theories. 
In this regard, of particular interest for us is the one with a non-minimally coupled scalar field because minimally 
coupled scalar source had already been investigated 0, 0|- Typical examples would be the Brans-Dicke theory @ and 
the the gravity with a dilaton field arising for instance in the string theory. 

In this paper, we extend the Horava-Lifshitz gravity to include the Brans-Dicke field as a concrete example of the 
non-minimally coupled scalar field. It turns out that such an extension is possible within the context of the detailed 
balance condition and in the IR limit reduces to the four-dimensional Brans-Dicke theory with negative cosmological 
constant when parameters of the theory satisfy certain conditions. 

We then study cosmological implication of the theory assuming homogeneity and isotropy. Without the negative 
cosmological constant and the dark radiation term[l(| the equations are those of the Brans-Dicke theory. So, we 
concentrate only on their effects on cosmology. Still, we find several interesting features. In the early universe limit 
there exists a solution where the scale factor a(t) grows like i 1 / 2 , which corresponds to the behavior of the universe 
in the presence of the normal radiation. Furthermore, in the large universe limit we find a solution which increases 
exponentially in spite of the existence of a negative cosmological term. This is contrary to the usual expectation that 
the exponential solution is possible only for a positive cosmological constant [111, EH- Both these aspects are possible 
because of the Brans-Dicke field. 

Let us consider the four-dimensional Brans-Dicke theory [8|, where the action is given by 

S = [ d^x^f~g (<pR - u^g^d^d^) . (1) 
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Decomposition of this action into (3 + 1) form, including the speed of light, c, yields 1 

V^9<PR * N^cj) (72 + c- 2 (K ab K ab - K 2 )) - 2N^c- 2 Kn - 2N^D 2 cj), (2) 
Vfyi<^0 = N^uj^c^tt 2 - N^Lo<j)- 1 D a (t)D a (j), (3) 

where the four metric g is decomposed into the lapse function N , the shift vector N a and the three metric q ab , and 
the corresponding three-dimensional covariant derivative and its scalar curvature are denoted respectively by D a , 72. 
The Brans-Dicke parameter is assumed positive, lo > 0. In the first equation an irrelevant total divergence term was 
dropped. The time derivatives of the three-metric and the scalar field are encoded in the following quantities; 

K ab = 7^{9ab- D a N b - D b N a ), (4) 

7T = l(^_tf°0 o 0). (5) 

Using the above result the Brans-Dicke action can be split into the two parts Sbd — S§ D + Sg D , where the kinetic 
and potential parts are 

S% D = c~ x J dtd 3 xN^/q (<j)(K ab K ab ~ K 2 ) - 2Ktt + w0~V) , (6) 

Sl D = c J dtd 3 xN^{<j)R-2D 2 (f)~LO(j)- 1 D a (j)D a (j)) . (7) 

Re-scaling the scalar field <f> and the corresponding field it, we find 

S§ D = J dtd 3 xN^((t){K ab K ab - K 2 )-2KTT + u(j)- 1 TT 2 ) , (8) 

S% D = c 2 J dt d 3 xNy^ (072 - 2D 2 (f> - uj(t)- 1 D a <j)D a (f>) . (9) 

Note that the factor of c 2 appears in front of the potential term. For the later purpose regarding the detailed balance 
it is important to express the kinetic part in the following matrix form; 



where 



" - / dt d* X NV- q ( K ab n ) ( -f[ ) ( K ; d ) , (10) 



Q abcd = ( q ac q bd + q ad q bc^ _ g ab q cd ^ 



Note that the matrix in the middle of the kinetic part of the action can be regarded as the supermetric on the space 
of (q a b, </>), naturally extending the DeWitt metric on the space of three-metrics. 

We intend to construct a Brans-Dicke type extension of Horava-Lifshitz gravity with the detailed balance condition. 
So, we choose the action of the form, Shlbd = Shlbd + ^hlbd-. where the kinetic part is 

SLbv = J dtfzNjq ( K ab tt ) ( ^J3 A ) -f\ ){ K ; d ) (12) 
and the potential part is of the form 



SZr.»n = - / dt*xNJq( % \% ) ^ d W J* [ £U ) ( 13 ) 
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for some suitable choice of function W(q,<j>). The supermetric G abcd {\) was slightly deformed compared to the Eq. 
(1111) to include the parameter A as usual, 

G abcd {\) = i (q ac q bd + q ad q bc ) - Xq ab q cd . (14) 



1 This result was considered in the context of conformal gravity in Ref. Il3ll . 
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Such a choice of the action is a natural generalization of the Horava-Lifshitz gravity in the context of the detailed 
balance condition. The factor of two was inserted in front of the variation of W with respect to <j> to compensate for 
different normalization in time derivatives in Eqs.(4) and (5). It is a straightforward matter to calculate the inverse 
supcrmctric. It comes out to be of form 



t> 1 Gabcd —Aq a b 

-Aq cd B</> 



(15) 



where 



with 



Qabcd = 2 (qa C qbd + qadqbc) - Mabqcd, (16) 



A - w (3A-l) + 3' (17) 
3A - 1 

B a,(3A-l) + 3' (18) 

— 1+wA . . 

= .(3A-D + 3 - (19) 

Note that this inverse supermetric is well-dehned even for A = 1/3 contrary to the pure gravity case and becomes 
singular instead when A = 1 and cj = — 3/2(oj>0is assumed in this work and it is nonsingular if A > 1/3.). The 
singualr case corresponds to the conformal scalar. If we take the limit of u> — > oo, A and B vanish and A = 1/(3A — 1), 
reproducing the pure gravity case. 
We choose 

W = ci J d 3 Xy^4>(TZ - 2A fc ) - c 2 J d i x^Lj(t>- 1 D a (j)D a (j). (20) 

In general all possible marginal and relevant terms can be included. The above choice of W corresponds to keeping 
only terms important in the infrared limit. Then, from 

— = - Cl A b H ab + Q ab , (21) 

oq a b 
1 SW 

2^ = ~ C ^ + Q > < 22 > 

where 

Q ab = ci (-4>G ab + D a D b cj) - q ab D 2 <p) , (23) 

Q = ci| - c 2 (-u^D 2 ^ + ^<t>- 2 D a <t>D a ct>) , (24) 

with G a b being the Einstein tensor constructed with the three-dimensional metric, we find after a straightforward 
calculation that 



Sl DHL = j dt d 3 xN^ + 0(4>K - ^-uj^D^Da^) + -f(-2D 2 cf>) 

- Jdt d 3 xN^d (Q ab cp- 1 g a bcdQ cd - 2AQ ab q ab Q + Bcj>Q 2 ) , (25) 



where 



2 3a; + 7 — 3A 



a = ^ ,(3A-l ) + 3 (26) 



cv + 5 — 3A 



' = -^ A ^(3A-l ) + 3 (2?) 



2( W + l)-a W (4-3A) 



? = ,(3A- l ) + 3 ■ <»> 
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Sbdhl\ir 

where 



The second line of Eq. (|25j) collects all the quadratic terms. 

When ci = c 2 and A = 1, the theory recovers four-dimensional difFeomorphism symmetry, as one can see from the 
fact that in the infrared limit the potential part of the action becomes 

-(ci) 2 A b ^±l ( dtd?xN^((} > (K-2k)-2D 2 cp-Lo ( j ) - 1 D a cj ) D a <j ) ) , (29) 

ZUI + 3 J 

3ui + 4 , , , 

= 25J+2) Ab ' (30) 

This expression coincides with that of the Brans-Dicke theory except that the cosmological constant term is present. 
Comparison with the kinetic part yields the speed of light 

As in the case of the Horava gravity the constant A& must be negative, consequently allowing only negative cosmological 
constant A. The Newton constant is related to the expectation value of the scalar field < > as follows, 

Gn = 1R (32) 
lD7r < <p > 

Now, we consider the homogeneous, isotropic cosmology. We will restrict ourself to the case of A = 1, c\ = C2, and 
set the speed of light to unity, i.e., c = 1. We choose vanishing shift vector iV a = 0, the three- metric to be the usual 
maximally symmetric ones with curvature constant k = —1, 0, +1, 

ds * = a 2 {t) ( + ^ 2 + ^2 2 A _ 

\ 1 — kr z I 



In this case the higher derivative terms become greatly simplified due to homogeneity and isotropy, 

1 

a* ~ a 

Substituting this result into the action, (I25p . yields the following mini-superspace action, 



ab = k Cl ^q ab , Q = 6k Cl ^. (34) 



Sbdhl = I dta 



_L L 6 ^f _ 6 ^ + u<jr x yf\ + N LB - 2A) + ^(fc Cl ) 2 4 
I\ \ a a I V a 2w + 3 a? 



(35) 



where the lapse field is set N = 1 after deriving the field equation. The field equations (including the matter) become 

,a. n a ■ 1 i-o ,3k , , 3w , 1 , 

3 - V + 3-0 - V + -J ~ A + fc Cl = -p ro , 36 

a a 2 ar 2(2w + 3) a 4 2 

-2-0 - - 2 - - 2-0 - - - A + fc Cl f\ = -p m , 37 

a a a 2 6(2w + 3) a 4 2 

w - i^-^ 2 + 3^0 - 3^0 - 3 A 2 - - A)0 - ^ (fc Cl ) 2 4 = °> ( 38 ) 
2 a a a a z 2(2u> + 3) a 4 

where the matter is assumed, for consistency, to satisfy the usual form of the continuity equation; 

Pm + 3-(pm+Pm) = 0- (39) 
a 

Only two equations are independent and can be chosen to be 

3H + 3H- --.(-) = -0 p „^ ( _- A) --(-), 

(2w + 3) (^+3*4) = i0- 1 ( /9ro -3p m )+2A+^, (40) 
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where H = (a /a) is the Hubble constant and 

R 2 3 " n. ,2 3^(3^ + 4) k 2 

B = 2Z7T3 (fcCl) = 2(c + 2) FA)' (41) 

The first in Eq. (1401) is the Friedmann equation of the Brans-Dicke theory with a negative cosmological term and 
the dark radiation term included. In the absence of those two terms the equations simply become those of the usual 
Brans-Dicke theory |l4|. Therefore, we restrict our attention to the new effects resulting from those two terms. 
Consider two limiting cases for vacuum, p m — p m = 0. First, for small a the dark radiation term dominates, so 

3** + 3*£-f„(i)° . (42) 



These equations can be solved by 



H = A, (44) 
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„2 ■ 



(45) 



with the two constants h and g satisfying 



3h 2 + 3hg-~ojg 2 = -\b 2 , (46) 

(2u + 3)g(g + h) = B 2 . (47) 

Eliminating the dark energy terms from the equations we find two possibilities h — —(\/2)g, or 3h = — (cj + 3)g. Only 
the first one gives rise to a solution when B 2 is positive, so we get 



h= --g = ±J-^— . (48) 

The two signs represent contracting and expanding phases. Solving for a and (f> reads 

a 2 (t) = 2ht, (49) 

m = ^, (50) 

where (f>o is the integration constant. This early universe behavior, a(t) ~ i 1 / 2 , is the one corresponding to the normal 
(not dark) radiation source. This is a rather unexpected result due to the scalar field. 

Secondly, in the large a limit cosmological term dominates over the curvature and dark radiation terms. We get 

SH' + SH^- 1 -^) 2 = A, (51) 

(2- + 3)((!)- + (|) 2 +3i/^ = 2A. (52) 

Similarly to the previous case, these two equations are solved by 

H = h, (53) 

\ - * ^ 

with h and g satisfying 

3h 2 + 3hg - ^ujg 2 = A, (55) 

(2w + 3)g(g + h) = A. (56) 



6 



In this case, we get either h = —\g or h — {lu + l)g. Only the first is allowed for negative A, and we find 



h = -i=±J^L. (57) 
2 V 2lu + 3 V ' 

For the positive sign the solution represents the universe exponentially expanding. It is interesting to note that such 
a solution exists even for a negative cosmological constant. This is in sharp contrast to the Horava gravity case. Note 
that the effective cosmological constant 3h 2 is suppressed by the factor of lu. 

To summarize, we have condtructcd a Brans-Dicke extension of the Horava-Lifshitz gravity with the detailed balance 
condition satisfied. We have investigated its IR limit and shown that the resulting IR theory is the Brans-Dicke theory 
with a negative cosmological constant and a dark radiation term. By studying its cosmological solutions we have shown 
that exponentially expanding solution at late time and power law expanding solution at early time can exist. This is 
in contrast with the pure gravity of Horava. Although we focused on the Brans-Dicke theory in this paper the analysis 
can be generalized to other non-minimally coupled scalar field gravity theory. It would be interesting to investigate 
further cosmological aspects of the resulting theories. 
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I. INTRODUCTION 



Recently, a new theory of gravity has been proposed by HoravapHIj]. This theory, being based on anisotropic 
scaling of space and time, breaks the spacetime symmetry. It has a much better UV behavior than the theories with 
the spacetime diffcomorphism symmetry, but expected to reduce to Einstein's gravity in the infrared limit recovering 
the spacetime diffeomorphism symmetry 1 . Physical constants such as the speed of light, Newton's constant, and 
cosmological constant all emerge from the relevant deformation of the non-relativistic theory at short distance. These 
interesting features as well as other related findings have received a great deal of attention^, @] ■ 

On the other hand, there are many alternative theories and extensions of the Einstein theory. In particular, various 
gravity models with scalar fields have been considered in the context of cosmology to explain the behaviors of the 
, universe in the early stage as well as in the late staged, 0. Therefore, it would be interesting to consider similar 
extensions in the context of Horava's theory. In this regard, of particular interest is the one with a non-minimally 
| coupled scalar field 2 , typical examples being the Brans-Dicke field[ll[ and the dilaton field [l2|. 
t-H . In this paper, we take the Brans-Dicke field as a concrete example of the non-minimally coupled scalar field and 
consider its inclusion into the framework of Horava's gravity. Originally, Horava introduced the concept of the detailed 
• *"H , balance condition as a way of reducing the choice of the potential, motivated by analogous methods used in quantum 
critical systems. Recently, many serious problems were reported|l3[ associated with strictly imposing this condition. 
However, some of the problems can be alleviated by softly breaking the condition. Although the fate of the detailed 
balance condition remains to be seen, it will be interesting to see if the detailed balance condition can be maintained 
when we try to non-minimally couple the scalar field to the Horava-Lifshitz gravity. 

It turns out that such an extension is possible and it reduces to the four-dimensional Brans-Dicke theory with 
negative cosmological constant when only the lowest order derivative terms are kept and parameters of the theory are 
chosen to satisfy certain conditions. 

We then study cosmological implication of the theory assuming homogeneity and isotropy, and including the 
curvature-squared terms. Because of the symmetries these higher order terms become a single term proportional 
to a~ 4 which can be regarded as the radiation with negative energy. However, it is not strictly so because the nor- 
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mal matter would couple with the inverse of the Brans-Dicke scalar field. We concentrate only on their effects on 
cosmology. We find several interesting features. We discuss them in Sec. 3. 



II. CONSTRUCTION OF THE MODEL 

Let us consider the four-dimensional Brans-Dicke theory where the action is given by 

S = J d^xj~g (<j)R - u^g^d^d^) . (1) 

Decomposition of this action into (3 + 1) form, including the speed of light, c, yields (See Ref. |l(|, for instance.) 

V^g<j)R ~ N^cj) (K + c- 2 (K ab K ab - K 2 )) - 2N^q~c7 2 K-K - 2N^D 2 <j), (2) 
-y/=gu)<jT x gi iV dp<l)d v <t> = Ny/qu^c^ir 2 - N^Lo<j)- 1 D a (t)D a (j), (3) 

where the four metric g is decomposed into the lapse function N, the shift vector N a and the three metric q ab , and 
the corresponding three-dimensional covariant derivative and its scalar curvature are denoted respectively by D a , 1Z. 
The Brans-Dicke parameter is assumed positive, uj > 0. In the first equation irrelevant total divergence terms were 
dropped out. The time derivatives of the three-metric and the scalar field are encoded in the following quantities; 

K ab = 7^(9ab- D a N b - D b N a ), (4) 

7T = ±(<j>-N a d a c(>). (5) 

Using the above result the Brans-Dicke action can be split into the two parts Sbd = &bd + ^BD' where the kinetic 
and potential parts can be written after re-scaling of the scalar field (j> and the corresponding field ir as 



S§ D = / dt (fxNy/d ((f>(K ab K ab - K 2 ) - 2K?r + lo^tt 2 ) , (6) 

S% D = c 2 J dt d 3 xNy^ (<j)K - 2D 2 cj) - uj(t)- 1 D a <j)D a (f>) . (7) 

Note the factor of c 2 in front of the potential term. For the later purpose regarding the detailed balance condition it 
is important to express the kinetic part in the following matrix form; 

S§ D =Jdt d?xN^~q (K ab n)( ~f[ ) ( K ; d ) , (8) 

where 

G abcd = ( q ac q bd + q ad q bc^ _ g ab q cd (g) 

The matrix in the middle of the kinetic part of the action can be regarded as the supermetric on the space of {q a bi 
naturally extending the DeWitt metric on the space of three-metrics. 

We intend to construct a Brans-Dicke type extension of Horava-Lifshitz gravity with the detailed balance condition. 
So, we choose the action of the form, Shlbd = 5f iBD + S^ LBD , where the kinetic part is 

S£ LBD = J dtd\N^ ( K a b 7T ) ( ^2 A) ) ( ^ ) (10) 

and the potential part is of the form 

SlLsn = - / dtfxNrf 0G :p A) -f[ ) 1 ( m ) (ID 

for some suitable choice of function W(q, tf>). The supermetric G abcd {\) was slightly deformed compared to the Eq. 
(|9]) to include the parameter A as usual, 

G ahcd {\) = - (q ac q bd + q ad q bc ) - Xq ab q cd . (12) 
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The factor of two was inserted in front of the variation of W with respect to <f> to compensate for different normalization 
in time derivatives in Eqs.(4) and (5). It is a straightforward matter to calculate the inverse supermetric. It comes 
out to be of form 



l Qabcd —Aq a b 
-Med B(j> 



(13) 



where 



with 



Gabcd — - (qacQbd + q a dqbc) - McMcd, (14) 



1 3A-1 - 1 + wA 

~ w(3A-l) + 3' w(3A-l) + 3' w(3A-l) + 3' (15) 

Note that this inverse supermetric is well-defined even for A = 1/3 contrary to the pure gravity case and becomes 
singular instead when A = (u> — 3)/3w, for instance when A = 1 and oj = —3/2 corresponding to the conformal scalar 
case (We assume w > in this work.). If we take the limit of lu — >• oo, A and B vanish and A = A/ (3A— 1), reproducing 
the pure gravity case. 
We choose 

W = ci J <P X yfq<j){n - 2A h ) - c 2 J d 3 x y fq~Lj(t)- 1 D a (f)D a (j). (16) 

In general all possible marginal and relevant terms can be included. The above choice of W corresponds to keeping 
only terms important in the infrared limit. Then, after a straightforward calculation Eq. (JTTJ) can be written as 

s hlbd = J dt d3xN V9 | a + P(4>K ~ ^-u)(j)~ X D a (j)D a (t>) + 7(-2L> 2 <; 

dt d 3 xN^d {Q ab cp- 1 g a bcdQ cd - 2AQ ab q ab Q + BtfiQ 2 ) , (17) 
where 

/ a \2 3w + 7 - 3A 
a = (GlAfe) o,(3A-l) + 3 (18) 

w(3A — 1) + 3 

, 2(cj + 1)-^w(4-3A) 

7 = -( C i) 2 A b ^ ' * \ '-, (20) 

uj(3\ — 1) + 3 

and 

Q afc = d ^-0(ft a6 - ^q Qb ) + £> a £> 6 (/> - q ah D 2 <y\ , (21) 
Q = (-uxT X D 2 (t> + ^cj>- 2 D a cj>D a c^ . (22) 

The second line of Eq. (jTT)) has quadratic terms only. 

When ci = c 2 and A = 1, the theory recovers four-dimensional diffeomorphism symmetry, as one can see from the 
fact that in the infrared limit the potential part of the action becomes 

SbdhlUr = ~( c i) 2Ab £^| j ' dtd 3 xN^q-(<b(K-2A)-2D 2 <j)-Lj^ 1 D a (bD a 4,) , (23) 



where 



= 2(a7T2) (24) 
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This expression coincides with that of the Brans-Dicke theory except that the cosmological constant term is present. 
Comparison with the kinetic part yields the speed of light 

As in the case of the Horava gravity the constant A& must be negative, consequently allowing only negative cosmological 
constant A. The Newton constant is related to the scalar field 4> as follows, 

GN = < 26 > 

III. COSMOLOGICAL SOLUTIONS 

Now, we consider the homogeneous, isotropic cosmology. We restrict ourself to the case of A = 1, c\ = C2, and set 
the speed of light to unity, i.e., c = 1. We choose vanishing shift vector N a = 0, and the three-metric to be the usual 
maximally symmetric ones with curvature constant k = —1, 0, +1, 

ds 2 = fl 2( t ) ( + r 2(^2 + gin 2 6d(f) 2-\ (1) 

\ 1 — kr 2 J 
In this case the higher derivative terms become greatly simplified due to homogeneity and isotropy, 

Q ab = kc x t i q a \ Q = 6k Cl \. (2) 
a 2 a 2 

The field equations become 

in* Mint 1 (K 2 3fc IB 2 

3H +3H-- -oj(-) = ~(j> R m ~ — +A--( — ), 



2 K (t>' T r a 2 2 a 

T 2 



<b u ,(/>,, o 1 1 k 1 , B 2 , 

-2H-3^-^-2 J ff^--(5) 2 = 0-V„ + --A--( — ), 

2 (x> 2 a z a* 



(2w + 3)^|+3ff|j = ^- 1 (p m -3 Pm ) + 2A+^, (3) 

together with the usual form of the continuity equation for the matter density p m for consistency, where H = (a /a) 
is the Hubble constant and 

2u + 3 y ' 2(w + 2) (-A) V ; 

The first equation in Eq. ([3]) is the Friedmann equation of the Brans-Dicke theory with a negative cosmological 
term and the dark radiation term included. Again, we emphasize that they do not have 0" 1 coupling in contrast to 
the normal matter. In the absence of those two terms the equations simply become those of the usual Brans-Dicke 
theory (T7). Therefore, we restrict our attention to the new effects resulting from those two terms. 

For simplicity, assume that the matter is absent, i.e., p m — p m — 0. First, consider the case where the dark 
radiation like term dominates, so Eq. ([3]) reduces to 



2 2 V a 4 

4> uj o 1, B 2 . 

-2H-3^-^-2F$--(|))» = --(^), 



ito + B) i + 3H!t) = 5. (5) 



a 4 



To further simplify these equation we set 



X = H + Y, Y=m), (6) 
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in terms of which, they can be written as 



where A = 2lu + 3. Setting 



Eq. © yields 



3X 2 -AY 2 
3 



X 



2a 4 ' 

^XY--Y 2 + -{^-), 
2 6 V 2a 4 ^' 

B 2 



AY = -3AXY + AY + — 

2 a 4 



(7) 



Vsx = 
Vay = 



2a 4 
B 2 
2a 4 




sinh^, 
cosh#, 



(8) 




, —cosh 9 H — sinh 

2a 4 Vv^ VA 



(9) 



Note that 9 is always negative for us > 0, which means that goes from negative infinity to positive infinity as time 
flows. Combine this result with 



to get 



= X-Y = 



dloga 
d6 



^^=cosh!# H =sinh 

y/A V3 



75 sinh 



-4=cosh # x 

VA 



-4= cosh 9 H — 7=sinh 

, V3 VA 



(10) 



(11) 



Although this equation can be integrated, the resulting expressions can be quite complicated, asymptotic behavior of 
a(t) at early and late times can be easily determined. A straightforward analysis shows that the scale factor vanishes 
at initial time and finite later time. The universe it describes expansion from a singularity and within a finite time it 
collapses. 

When the cosmological term dominates, Eq. becomes 



3A 2 -AY 2 
3 



A. 



X 



2 X > 



A , 

XY Y 2 

2 



AY = -3AXY + AY 2 



A. 



A 
2 



(12) 



With 



V3X = V^A sinh 9, 
y/AY = V^A cosh 9, 



(13) 



we find a slightly different equation for 9, 



-A f\/3cosh 9 - -F| smh 9 



(14) 



General behavior of the solution is the same as the previous case. 
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IV. CONCLUSION AND DISCUSSION 

To summarize, we constructed a Brans-Dicke type extension of the Horava-Lifshitz gravity maintaining the detailed 
balance condition. Although strict imposition of the detailed balance condition is known to have many problems, 
one can either break the condition or simply treat the resulting terms as an important contribution to the potential. 
We have not discussed the issue of projectability in this work. At this level, our model can be incorporated into any 
version. 

Furthermore, We investigated its low energy limit and shown that the resulting theory is the Brans-Dicke theory 
with a negative cosmological constant. In Brans-Dicke theory one can incorporate cosmological constant term in two 
ways. One is treating it as a vacuum expectation value from the matter sector and the other is what we have done in 
this work. 

We studied the curvature contribution up to quadratic order in the context of homogeneous and isotropic cosmology. 
The resulting theory is the Brans-Dicke theory with a negative cosmological constant and a radiation-like term with 
a negative energy density. They are somewhat different from usual matter in that they do not have coupling. We 
analyzed their effects and showed that the resulting solution has the general behavior of big rip. This is in contrast 
with the pure gravity case of Horava. 

Although we focused on the Brans-Dicke theory in this paper the analysis can be generalized to other non-minimally 
coupled scalar field gravity theory. It would be interesting to further investigate cosmological aspects of the resulting 
theories. 



Acknowledgments 

THL was supported by the Soongsil University Research Fund. PO was supported by the National Research 
Foundation of Korea(NRF) grant funded by the Korea government (MEST) through the Center for Quantum Space- 
time(CQUeST) of Sogang University with grant number 2005-0049409. 



[1] P. Horava, JEEP 0903 (2009) 020. 

[2] P. Horava, Phys. Rev. D79 (2009) 084008. 

[3] P. Horava, Phys. Rev. Lett. 102 (2009) 161301. 

[4] D. Bias, O. Pujolas and S. Si biryakov, Phys. R ev. Lett. 104 (2010) 181302, Phys. Lett. B688 (2010) 350. 
[5] J. Bellorin and A. Restuccia, larXiv:1004.0055l 
[6] G. Calcagni, JEEP 09, (2009) 112. 

[7] H. Lu, J. Mei, and C. N. Pope, Phys. Rev. Lett. 103, (2009) 091301; E. Kiritsis and G. Kofinas, Nucl. Phys. B821, (2009) 

467; T. Takahashi and J. Soda, Phys. Rev. Lett. 102, (2009) 231301; R. Brandenberger, Phys. Rev. D80, (2009) 043516; C. 

Charmousis, G. Niz, A. Padilla, and P. M. Saffin, JEEP 0908, (2009) 070; S. Mukohyama, JCAP 0906, (2009) 001; R.-G. 

Cai, L.-M. Cao, and N. Ohta, Phys. Rev. D80, (2009) 024003; H. Nastase, larXiv:0904.3604l A. Kehagias and K. Sfetsos, 

Phys. Lett. B678, (2009) 123; X. Gao, arXiv:0904.4"T87l M. Li and Y. Pang, JEEP 0908, (2009) 015; T. P. Sotiriou, M. 

Visser, and S. Weinfurtner, JEEP 0910, (2009) 033; R.-G. Cai, Y. Liu, and Y.-W. Sun, JEEP 0906, (2009) 010; Y.-S. 

Piao, Phys. Lett. B681, (2009) 1; R.-G. Cai, B. Hu, and H.-B. Zhang, Phys. Rev. D80, (2009) 041501; R.-G. Cai, L.-M. 

Cao, and N. Oht a, Phys. Lett. B6 79, (2009) 504; E. O Colgain and H. Yavartanoo, JEEP 0908, (2009) 021; Y. S. Myung 

and Y.-W. Kim, larXiv:0905.0T79l M. Visser, Phy s. Rev. D 80, (2009) 025011; D. Orlando and S. Relfert, Class. Quant. 

Grav. 26, (2009) 155021; A. Ghodsi and E. Hatefi. larXiv:0906. 12371 R. A. Konoplya, Phys. Lett. B679, (2009) 499; B. R. 

Majhi, Phys. Lett. B686, (201 0) 49; M. R. Setare and M. Jamil, JCAP 1002, (2010) 010; M. Jamil, E. N. Saridakis and 

M. R. Setare. larXiv: 1003.08761 Y. Cai and E. N. Saridakis, JCAP 0910, (2009) 020; S. Dutta and E. N. Saridakis, JCAP 

1001, (2010) 013; M. I. Park, JEEP 0909, (2009) 123. 
[8] B. Ratra and P.J.E. Peebles, Phys. Rev. D37 (1988) 3406; RR Caldwell, R. Dave and P.J. Steinhardt, Phys. Rev. Lett. 

80, (1998) 1582. 

[9] C. Armendariz-Picon, V. Mukhanov and P.J. Steinhardt, Phys. Rev. Lett. 85, (2000) 4438; T. Okabe and M. Yamaguchi, 
Phys. Rev. D62, (2000) 023511; M. Malquarti, E.J. Copland, AR. Liddle and M. Trodden, Phys. Rev. D67, (2003) 123503. 
[10] S. Mukohyama, JCAP 0906, (2009) 001; E. Kiritsis and G. Kofinas, Nucl. Phys. B821, (2009) 467; B. Chen and Qing-Guo 

Huang, Phys. Lett. B683, (2010) 108. 
[11] C. Brans and R. H. Dicke, Phys. Rev. 124, (1961) 925. 

[12] C.G. Callan, D. Friedan, E.J. Martinec and M.J. Perry, Nucl. Phys. B262, (1985) 597. 

[13] G. Calcagni, Phys. Rev. D81 (2010) 044006; C. Bogdanos and E. N. Saridakis, Class. Quant. Grav. 27 (2010) 075005; T. 
Sotiriou, M. Visser and S. Weinfurtner, Phys. Rev. Lett. 102, (2009) 251601. 

[14] K. Uehara, CW. Kim, Phys. Rev. D26, (1982) 2575. 

[15] M. Arik, M. C. Calik, JCAP 0501, (2005) 013, |gr-qc/0403108 



7 



[16] B. Kelleher, Class. Quantum Grav. 21, (2004) 483. 

[17] Hongsu Kim, Mon. Not. Roy. Astron. Soc. 364, (2005) 813. 



